TABLE  OF  CONTENTS 


Page 

ABSTRACT  .  1 

I.  INTRODUCTION  .  2 

II.  ACCOMPLISHMENTS . 

A.  Papers  published .  ** 

B.  Fapers  accepted  for  publication  .  5 

C.  Papers  submitted  .  5 

III.  MEASUREMENTS  AND  RESULTS  .  6 

A.  Specimens .  6 

B.  Pressure  and  Temperature  'easurement  .  7 

C.  Velocity  Measurements  ar  Elastic 

Properties  Determinatior  .  ° 

D.  Experimental  Results  .  H 

IV.  RECOMMENDATION  FOR  FUTURE  WORK . 18 

References . 2i 

APPENDIX  A 

ABSTRACTS  OF  PAPERS . 23 

APPENDIX  B 

MEASUREMENT  OF  SOUND  VELOCITY  BY  MEANS  OF 
ULTRASONIC  INTERFEROMETRY  .  31 

A.  Phase  Comparison  .  31 

B.  Pulse  Superposition  .  35 

References . 39 

APPENDIX  C 

ADIABATIC  TO  ISOTHERMAL  CORRECTIONS  TO  THE 

SECOND  ORDER  .  40 

Reference  .  ^ 


TABLE  OF  CONTENTS  (continued) 


Page 

APPENDIX  D 

EXPERIMENTAL  DATA . 45 

A.  MgO  .  . . 45 

B.  A1203  . 

References  .....  .  47 


LIST  OP  TABLES 


following 

page 


I.  Compression  of  properties  determined  at 

-78 . 5°C  and  3^ .  6°C  for  MgO .  12 

II.  Comparison  of  properties  determined  at 

-78.5°C  and  25°C  for  A1203  .  12 

III.  Comparison  of  calculated  with  measured 

compression  for  MgO .  15 

IV.  Comparison  of  calculated  with  measured 

compression  for  A1203  15 

APPENDIX  C 

I.  Correction  terms  in  the  calculation  of  (3Brp/3P),p  43 

II.  Corrections  for  computing  (3Bm/3P)  of  Mgo 

at  -78  5°C . T .  ^3 

III.  Corrections  for  computing  (3Sm/3P)  of 

A1203  at  25°C . T .  43 

IV.  Corrections  for  computing  (3Bm/3P)  of 

A1203  at  -78 . 5°C . ! .  43 

APPENDIX  D 

I.  Longitudinal  and  shear  velocities  with 

temperature  at  1  atm  for  MgO .  45 

II.  Longitudinal  and  shear  velocities  with 

pressure  at  -78.5°C  and  34.6°C  .  46 

III.  Variation  of  adiabatic  bulk  modulus  with 

pressure  at  34.6°C  and  -78.5°C  for  MgO .  46 

IV.  Longitudinal  and  shear  velocity  as  function 

of  temperature  at  1  atm .  46 

V.  Longitudinal  and  shear  velocity  as  function 

of  pressure  at  25°C  and  -78.5°C .  46 

VI.  Variation  of  adiabatic  bulk  modulus  with 

pressure  at  25°C  and  -78.5°C  for  A12C>3 .  46 


LIST  OF  FIGURES 

following 

page 


1.  Block  Diagram  of  the  4  Kilobar  Pressure 

System . 7 

2A.  T-P  Plot  at  Constant  Sound  Velocity  for 

MgO  with  Geothermic  Gradients  Superposed  ....  17 

2B.  T-P  Plot  at  Constant  Sound  Velocity  for 

with  Geothermic  Gradients  Superposed  ...  17 

3.  Velocity-Depth  Profiles  Obtained  for  MgO 

and  A120^  .  . .  17 

APPENDIX  B 

IA.  Ultrasonic  Interferometer  for  Phase 

Comparison  .  . .  31 

IB.  Ultrasonic  Interferometer  for  Pulse 

Superposition  .  35 

2,  Illustrating  the  Meaning  of  the  Integer  N  .  .  .  37 


-  1  - 

ABSTRACT.  The  techniques  of  ultrasonic  interferometry  were 
used  to  measure  the  isotropic  sound  velocities  and  their 
derivatives  for  polycrystalline  specimens  of  aluminum  oxide 
and  magnesium  oxide.  The  pressure  derivatives  were  deter¬ 
mined  to  4  kbars  at  room  temperature  and  at  the  -78.5°C. 

The  temperature  derivatives  were  measured  at  1  atmosphere. 
The  velocities  and  their  derivatives  at  room  temperature 
are 


MgO 

A  1^0 

VP 

9-7711 

10.845 

km/sec 

vs 

5.9682 

6.3730 

km/sec 

dvp/dP 

7 . 711*10“3 

5 . 175*10“3 

km/sec/kbar 

dv3/dP 

4.351*10“- 

2 . 207*10- 3 

km/sec/kbar 

dVp/dT 

-5.0  xicr4 

-3-7  *10“4 

km/sec/°C 

dvs/dT 

-4.8  *10“4 

-3.1  *10"4 

km/sec/°C 

The 

isothermal  pressure 

derivatives  of 

the  bulk  modulus 

were  determined  to  be  (dBT/dP)^  Q  =  3-94  and  (dB(p/dP)A1  0  = 


2  3 


3.99  at  room  temperature. 

Prom  these  data,  the  critical  temperature  gradient  for 
velocities,  (dT/dP)v,  were  evaluated,  and  the  velocity  behav¬ 
ior  of  these  materials  as  a  function  of  temperature  and 
depth  in  the  mantle  are  discussed.  It  was  found  that,  because 

(dT/dP)v  >  (dT/dP)v  ,  they  exhibit  a  shear  velocity  minimum 
VP  vs 

and  a  less  pronounced  longitudinal  velocity  minimum  under 
conditions  likely  to  exist  in  the  upper  mantle. 


la 


An  expression  was  derived,  employing  the  measured 
values  of  the  hulk  modulus  B0  and  its  pressure  derivative 
B0',  which  predicts  compression  at  extremely  high  pressures. 
Using  the  derived  expression 

ln(V0/V)  =  (1/Bq’)  in  CBq‘ (P/Bq)  +  i] 

and  the  vr.lues  of  B0  and  BQ'  for  MgO  and  tl2°3  reported  here, 
the  agreement  between  the  calculated  compression  and  measured 
compression  reported  in  the  literature  was  within  0.5$  to 
350  kilobars. 
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I .  INTRODUCTION 

The  objective  of  this  research  program  was  to  undertake 
investigation  of  the  behavior  of  the  pressure  and  tempera¬ 
ture  derivatives  of  homogeneous  mineral  specimens;  the  data 
obtained  are  applicable  to  the  resolution  of  questions 
regarding  the  behavior  of  seismic  velocities  in  a  homoge¬ 
neous  mantle.  Some  of  the  problems  to  which  results  of 
these  studies  apply  are  the  lateral  variations  of  velocity 
at  constant  depth;  variation  of  seismic  velocities  with  depth; 
dispersion  of  seismic  velocities  with  temperature;  and, 
finally,  the  problem  of  the  'shadow  zone'  phenomenon  in  the 
upper  mantle. 

The  program  differs  from  the  usual  experiments  on  the 
sound  velocity  in  rocks  in  several  ways.  First,  the  mate¬ 
rials  employed  were  synthetically  prepared  polycrystalline 
mineral  specimens.  The  specimens  were  prepared  in  a  manner 
so  that  they  were  of  essentially  zero  porosity  and  controlled 
(small)  grain  size.  Second,  the  technique  used  for  the  mea¬ 
surement  of  the  sound  velocities  and  their  derivatives  was 
ultrasonic  interferometry,  which  permitted  accurate  deter¬ 
mination  of  the  velocity  derivatives  in  the  P-T  plane.  This 
accuracy  was  sufficient  to  permit  the  calculation  of  the 
compression  of  the  material  well  beyond  the  range  of  the 
actual  measurements. 
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Investigations  on  the  behavior  of  rock  specimens  have 
taken  two  approaches.  These  are  the  direct  measurement  of 
velocity  with  pressure  or  temperature ; 1 » 2 > 3  or  the  direct 
measure  of  the  compressibility.4  In  both  instances  the 
presence  of  pores,  cracks,  and  mineral  orientations  compli¬ 
cate  the  results.  Typically,  the  initial  slopes  represent 
extrinsic  rather  than  intrinsic  effects;  the  1  atm  veloc¬ 
ities  or  compressibilities  are  obtained  by  extrapolating 
back  from  higher  pressure  data  where  the  intrinsic  behavior 
is  observed.  Further,  the  techniques  are  generally  not 
sufficiently  accurate  to  define  the  pressure  derivatives. 
Temperature  measurements  on  rocks  suffer  the  further  disad¬ 
vantage  of  the  irreversible  chemical  changes  within  the  rock 
when  It  is  heated,  so  that  low  temperature  measurements 
would  not  correspond  to  high  temperature  measurements  on 
the  same  materials.  Measurements  on  mono-mineral  specimens 
of  extremely  good  quality  do  not  suffer  these  disadvantages. 
In  particular,  intrinsic  properties  are  measured  over  the 
entire  experimental  range,  and  with  sufficient  accuracy  to 
define  the  derivatives.  Thus,  the  results  of  these  experi¬ 
ments  leads  to  precise  numerical  data  applicable  to  the  pro¬ 
pagation  of  sound  velocity  and  the  elastic  properties  within 


the  earth. 
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II.  ACCOMPLISHMENTS 

We  have  measured  the  velocities  and  their  pressure  and 
temperature  derivatives  of  two  polycrystalline  specimens, 
one  of  MgO  and  the  other  A^O^.  The  results  have  been  pub¬ 
lished  in  a  series  of  papers  and  their  applicability  to  the 
problems  of  the  deep  earth  discussed.  Analytic  methods  for 
computing  the  compression  and  velocities  at  extremely  high 
pressures,  employing  the  data  obtained  at  relatively  low 
pressure  were  developed.  In  all,  ten  papers  were  produced. 
Seven  have  been  published,  two  have  been  accepted  and  are 
in  press,  and  one  has  been  submitted.  The  titles  of  tnest 
papers  are  listed  below.  Abstracts  of  these  papers  appear 
in  Appendix  A. 

A .  Papers  published: 

1.  "Conditions  for  a  Density  Minimum  in  the  Upper 
Mantle,"  0.  L.  Anderson,  J .  Geophys .  Res . ,  70[6],  1457-1461, 
1965. 

2.  "The  Relation  between  Refractive  Index  and  Density 
of  Minerals  Related  to  the  Earth's  Mantle,"  0.  L.  Anderson 
and  E.  Schreiber,  J .  Geophys .  Res .  ,  £0[6],  1463-1471,  1965. 

3.  "Lattice  Dynamics  in  Geophysics,"  0.  L.  Anderson, 
Trans.  N.Y.  Acad.  Sol.  .  Ser.  II,  2J_C 33 ,  298-303,  1965. 

4.  "The  Bulk  Modulus-Volume  Relationship  for  Oxide 
Compounds  and  Related  Geophysical  Problems,"  0.  L.  / iderson 
and  J.  E.  Nafe,  J .  Geophys ,  Res .  ,  7C)[l6],  3951—3963 >  1965. 
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5.  !lAn  Approximate  Method  of  Estimating  Shear  Velocity 
from  Specific  Heat,"  0.  L .  Anderson,  (Letter)  J ■  Geophys .  = 

Res^,  jot  10],  14726-4728 ,  1965. 

6.  "The  Pressure  Derivative  of  the  Sound  Velocities  of 
Polycrystalline  Magnesia,"  0.  L.  Anderson  and  E.  Schreiber, 

J.  Geophys.  Fes.,  70[20],  5241-5248 ,  1965. 

7.  "Two  Methods  for  Estimating  Compression  and  Sound 
Velocity  at  Very  High  Pressures,"  Proc.  Nat.  Acad.  Scl., 
5iC3],  667-673,  1965. 

B.  Papers  accepted  for  publication: 

1.  "The  Pressure  Derivatives  of  the  Sound  Velocities 

of  Polycrystalline  Alumina,"  E.  Schreiber  and  0.  L.  Anderson, 
J.  Am.  Ceram.  Soc . ,  in  press. 

2.  "Seismic  Parameter  <j»:  Computation  at  Very  High 
Pressure  from  Laboratory  Data,"  0.  L.  Anderson,  Eull.  Seism. 
Soc .  Am. ,  in  press. 

C •  Papers  submitted: 

1.  "Temperature  Dependence  of  the  Velocity  Derivatives 
of  Periclase,"  E.  Schreiber  and  0.  L.  Anderson. 
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III.  MEASUREMENTS  AND  RESULTS 


A  specimen  of  hot-pressed  magnesium  oxide  was  obtained 
from  the  ceramics  laboratory  at  the  University  of  California 
at  Berkeley.  This  material  was  99*9%  theoretical  density. 

The  aluminum  oxide  (lucalox)  was  obtained  from  the  General 
Electric  Company.  This  specimen  had  a  density  that  was 
99-15%  of  theoretical.  The  specimens  were  prepared  by  having 
two  basal  surfaces  ground  and  polished  parallel  to  within 
0.0001"  per  inch  and  flat  to  within  0.1  wavelength.  Specimen 
length  was  measured  with  a  micrometer  which  read  directly  to 
0.0001  cm.  The  micrometer  and  specimen  were  thermostated  on 
a  large  bloc'.c  of  aluminum.  The  micrometer  was  checked 
against  grade  A-l  gauge  blocks. 

Specimen  density  was  determined  using  the  method  of 
Archimedes.  This  was  chosen  a*  the  means  for  obtaining  the 
most  accurate  density  values.  This  is  important,  as  the 
elascic  properties  are  calculated  from  the  sound  velocities; 
and,  to  minimize  degrading  the  data  computed  from  the  sound 
velocities,  the  density  must  be  known  with  equivalent  accu¬ 
racy.  To  obtain  the  desired  accuracy,  it  is  necessary  to 
correct  for  all  bouyancy  effects  and  for  surface  tension 
forces  on  the  suspending  fiber.  To  minimize  the  latter 
effect,  it  is  desirable  to  use  a  fiber  which  has  a  large 
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contact  angle  at  the  air-water  interface,  and  as  small  a 
cross  section  as  possible.  Nylon  fiber,  4  mil  diamter,  and 
a  contact  angle  of  70°?0'  were  used.  A  Mettler  semi-micro 
balance,  reproducible  to  0.00003  gms,  was  employed  in  all 
the  mass  determinations.  It  should  be  noted  that  accurate 
bulk  densities  can  be  obtained  in  this  manner  only  for 
impermeable  specimens. 

B.  Pressure  and  Temperature  Measurement 

The  4  kbar  pressure  system  is  shown  in  Figure  1.  A 
60,000  psi  S-C  pump  is  the  primary  pressure  generator.  It 
delivers  oil  to  one  side  of  the  mercury  piston  gas  compressor. 
The  mercury  column  functions  as  a  frictionless  piston  and 
serves  to  separate  the  oil  and  gas  sides  of  the  system. 
Pressure  is  measured  on  a  0-75,000  psi  Heise  gauge,  which 
is  connected  to  a  dead  weight  tester  for  calibration.  The 
compressed  gas  is  fed  into  the  vessel  which  contains  the 
specimen  and  holder.  Leads  through  the  cap  of  the  vessel 
provide  the  means  of  applying  the  r-f  signal.  A  copper- 
constantan  thermocouple  is  brought  through  the  cap  of  the 
vessel  also.  In  this  manner,  the  actual  specimen  tempera¬ 
ture  is  monitored  within  the  pressure  vessel.  A  Leeds  & 
Northrup  K-3  potentiometer  was  used  to  measure  thermocouple 
output . 

For  the  low  temperature  runs,  the  vessel  was  placed 
in  a  cryostat  containing  crushed  dry  ice.  For  room 
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temperature  runs,  the  vessel  was  circled  with  copper  tubing 
through  which  water  was  circulated  from  a  thermostated  bath. 
Measurements  at  intermediate  temperatures  were  obtained  as 
the  vessel  was  slowly  cooled  to  the  temperature  of  dry  ice, 

C '  ¥.el-oclty  Measurements  and  Elastic  Properties  Determination 

The  techniques  for  measuring  the  sound  velocity  v/ere 
ultrasonic  interferometry.  The  methods  of  phase  comparison 
and  pulse  superposition  were  both  employed.  The  details  of 
these  methods  are  discussed  in  Appendix  B.  Both  methods 
depend  on  interfering  the  radio  frequency  energy  in  a  pulse 
applied  to  the  specimen  with  the  energy  in  the  return  echo. 

By  taking  advantage  of  the  interference  effects,  measurements 
of  transit  time  accurate  to  0.001$  are  attainable,  per¬ 
mitting  velocity  determinations  accurately  to  0.01$.  The 
sensitivity  of  the  method  for  detecting  changes  in  transit 
time  as  small  as  0.1  nanosecond  make  the  methods  ideal  for 
determining  velocity  change  as  a  function  of  pressure  or 
temperature.  In  practice,  one  measures  the  initial  velocity, 
v0j  at  ambient  conditions,  and  then  measures  the  change  in 
frequency,  (time"1 ) ,  as  a  function  of  pressure  or  temperature. 
The  velocity  at  the  new  condition  is  calculated  from 


V  =  V°  I 

\  °/  \  o  J 


(1) 
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where  f/f  is  the  ratio  of  frequency  at  the  new  condition 
to  the  frequency  at  the  initial  (ambient)  conditions  and 
i/i0  is  the  length  ratio  for  the  two  conditions  in  the  direc¬ 
tion  of  sound  propagation.  In  temperature  runs,  the  change 
in  length  is  determined  from  thermal  expansion  data.  For 
pressure  runs,  the  ratio,  £/£0,  is  obtained  by  integrating 
the  frequency  dependence  of  the  velocity  modes  with  pressure 
(see  Appendix  B).  In  this  manner  the  variation  of  the 
velocities  with  pressure  or  temperature  are  accurately 
determined.  Because  one  measures  the  velocities  directly 
and  the  density  is  determinable  [p  =  p0(A0/A)3],  the  adia¬ 
batic  elastic  properties  are  defined.  Of  particular  interest 
are  the  shear  modulus,  Gg,  bulk  modulus,  Ba,  and  Poisson's 
ratio,  a,  which  are  given  by 


a  pvs2 

(2a.) 

"  p(vp2  -  V3  Vs2) 

(2b) 

-  1/2  (i  -  c(vp/vs)2  -  lr1} 

(2c) 

where  vs  and  vp  are  the  shear  and  compression  velocities 
of  an  isotropic  solid. 

The  variation  of  the  moduli  with  pressure  or  temperature 
are  given  by 

(  3GS\ 

\a  T/p  "  20pVs<9vs/9T>P  “  Gs°V 


(3a) 


10 


20p 


3  vs 


+  1  +  a  yT 
v 


In  the  above  equations  p  is  in  gms/cm3,  ay  Is  the  volume 
expansion,  and  y  Is  the  Gruneisen  constant  given  by 


v  Bs 

7? 


(5) 


where  Cp  Is  the  heat  capacity  at  constant  pressure. 

To  compare  results  of  compressibility  obtained  from 
acoustic  data  with  those  obtained  from  compression  measure¬ 
ments,  one  must  take  Into  account  the  adiabatic  to  isothermal 
correction.  As  Gg  =  G^,  (the  T  subscript  is  used  to  represent 
the  isothermal  condition),  corrections  are  needed  only  for 
the  bulk  modulus.  This  may  be  obtained  through  the  relations 
(see  Appendix  C  for  details): 


Bt  -  3S(1  +  Tavy) 


(6) 
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The  isothermal  pressure  derivative  of  the  bulk  modulus 
is  given  by 


and  the  temperature  derivative  of  the  isothermal  bulk  modulus 
obtained  from  equation  (6)  is  given  below. 


D.  Experimental  Results 

Velocity  measurements  were  performed  on  the  MgO  and 
Al^O^  specimens  as  a  function  of  pressure  and  of  temperature, 
as  described  previously  and  in  Appendix  A.  These  measure¬ 
ments  were  performed  under  the  following  conditions: 

(1)  isothermal  pressure  runs  at  ambient  or  near  ambient 
temperature;  (k)  isobaric  temperature  runs  down  to  -78.5°C 
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(the  sublimation  temperature  of  dry  ice)  at  1  atmosphere 
pressure;  and  (3)  isothermal  pressure  runs  at  -78.5°C.  The 
detailed  results,  which  we  submitted  for  publication, 6>7»8 
are  given  in  Appendix  D.  The  pressure  runs  were  to  4  kbars, 
except  for  MgO  at  -78. 5°C .which  was  to  2  kbars.  In  Tables 
I  and  II  the  properties  computed  from  these  basic  data  are 
listed  for  magnesium  oxide  and  aluminum  oxide  respectively. 

The  data  indicates  that  the  shear  and  longitudinal 
velocities  of  both  materials  are  linear  with  pressure;  and  the 
bulk  modulus  is  linear  with  pressure  for  both  these  materials. 
This  is  important  because  we  are  now  able  to  use  two  methods9 
for  calculating  compression  into  the  hundred  kilobar  region 
from  this  precise  data  which  was  measured  at  low  temperatures . 

The  first  method  for  deriving  a  compression  curve 
depends  upon  the  approximation  that  the  isothermal  bulk 
modulus  is  linear  with  pressure;  i.  e., 


B  =  -V 


§  h*o  +  v  p 


(9) 


Expanding  the  volume  isothermally  by  a  Maclaurin's 
series  to  order  P4,  we  have 


V(P)  -  V(0)  +  V'(0)P  +|v»  1  (0)P2  +|v»'(0)P3  +^jrV""(0)P4  (10) 


•V 


TABLE  I.  Comparison  of  properties  determined  at  -78.5°C  and 


34.6°C  for  MgO 


PROPERTY 

TEMPERATURE 

UNITS 

-78.5°C 

34 ,6°C 

Long.  Velocity, 

Vp  at  P=0 

9.8179 

5,7662 

km/sec 

Shear  Velocity, 
vs  at  P=0 

6.0136 

5.9635 

km/sec 

(  3  Vp/ 8  P  )ip 

7.57*10’3 

7 .7HX10”3 

km/sec/kbar 

(3vs/3P)t 

4.07x10“3 

4.351xio”3 

km/sec/kbar 

( ^Vp/3T)p_0 

4.3X10"4 

5.0*10~4 

km/sec/°C 

(3vg/3T)p=Q 

3 .  6xl0~4 

4.8xlo-4 

km/sec/°C 

Adiabatic  Bulk  Mod. 
Bs,  P=0 

1729.4 

1717.0 

kbars 

Isothermal  Bulk  Mod. 

V  ?«o 

1717.1 

1691.0 

kbars 

(3B5/8P)t 

4.00 

3.92 

( 3Bip/3P 

4.02 

3.94 

( 3Bg/3T )p-Q 

-0.11 

-0.13 

kbars/°C 

(3BT/3T)p.0 

-0.25 

-0.29 

kbars/°C 

Gruneisen  Constant 

1.67 

1.60 

Poisson's  Ratio 

0.200 

0.203 

Density 

3.5903 

3.5800 

gm/CC 

TABLE  II.  Comparison  of  properties  determined  at  -78.5°C 


and  25°C  for  A1203. 


PROPERTY 

TEMPERATURE 

UNITS 

-78.5°C 

?5°C 

Long.  Velocity, 

V  P*0 

10.880 

10.845 

km/sec 

Shear  Velocity 
va  at  P=0 

w 

6.4028 

6.3730 

km/sec 

(3Vp/3P)T 

5 .075x10“3 

5 . 175x10~3 

km/sec/kbar 

Ov3/3P)t 

2.1O7><10“3 

2 . 207x10“3 

km/sec/kbar 

<3V3T)p*o 

-2 . 92*10“4 

-4 . 45*io~4 

km/sec/°C 

(3V8T)p.o 

-2 .07x10“4 

-3 .19x10“4 

km/sec/°C 

Adiabatic  Bulk  Mod. 

V  p=0 

2534.2 

2520.6 

kbars 

Isothermal  Bulk  Mod. 
Brp,  P=0 

2526.8 

2504.5 

kbars 

OB8/3P)t 

4.00 

3-98 

OBt/3P)t 

4.01 

3.99 

(SBs/3T)p 

-0.118 

-0.138 

kbar/°C 

(3BT/3T)p 

-0.163 

-0.207 

kbar/°C 

Gruneisen  Constant 

1.23 

1.32 

Poisson's  Ratio 

0.2351 

0.2363 

Density 

3.9773 

3.9720 

gms/cm5 
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where  the  primes  denote  derivatives  with  respect  to  pressure. 
Solving  for  the  compression  equation  using  (9)  ar.d  the  neces¬ 
sary  operations  on  (2a),  we  have 


where 

m  =  1/2(1  +  B  * ) 

n  =  1/6  [1  +  3B0»  +  2(B0')2  -  B0B0"] 
q  =  1/24  [ll(Bo')2  +  6Bq  '  +  6(Bq')H  1 

~  ^oBo”  “  6BoBo'Bo"  +  W"3 


(U) 


(12) 

(13) 


(HO 


Equation  (11)  holds  to  a  pressure  given  by 

I  *•' 

P*  «  3/2[B/(l  +2B0’)]  05) 


It  turns  out  that  P*  is  290  kbars  for  MgO  and  419  kbars  for 
AI2O2 • 

In  (13)  and  (14)  all  terms  Involving  B  * ’  and  B  * 1 '  are 

0  o 

ignored  consistent  with  (9)«  This  is  an  arbitrary  assumption 
which  can  be  refined  by  future  theoretical  and  experimental 
work. 

The  second  method  for  deriving  a  compression  curve  may 
be  found  by  assuming  that  the  Instantaneous  bulk  modulus  is 
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a  linear  function  of  pressure.  That  iB,  Instead  of  (9), 
we  have 


B 


b  +  aP 


(15) 


Murnaghan10  defined  ,nd  Integrated  (15)  to  find 


P  = 


b 

a 


(16) 


Equation  (16)  is  called  the  Murnaghan  logarithmic  equation. 
Murnaghan10  found  that  the  acoustic  value  of  a  is  B  ’•  and 
iu  appears  obvious  that  b  =  BQ.  Thus,  acoustic  data  defines 
the  parameters  used  in  the  Murnaghan  logarithmic  compression 
equation 


Using  the  isothermal  value;  ;f  the  bulk  modulus  and  its 
pressure  derivative,  we  have,  for  MgO  and  A12C>3,  the 
equations 


B0  =  1692  kbars 


(18a) 
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To  test  these  equation!,  values  of  (V/V  )  were  computed 
from,  equations  (18)  and  (19).  The  results  for  MgO  are  com¬ 
pared  with  the  compression  data  of  Perez-Albuerne  and 
Drickamer1 1  in  Table  III.  The  result  of  the  calculations 
for  A120^  are  compared  with  the  compression  data  of  Hart 
and  Drickamer12  and  are  shown  in  Table  IV. 

Another  Important  result  is  the  evaluation  of  the 
critical  temperature  gradient  for  MgO  and  A120^.  The  critical 
temperature  gradient  is  that  value  of  (dT/dP)v  which  when 
applied  to  the  material  just  causes  the  sound  velocity  to 
be  refracted  with  positive  curvature.  Thus,  If  the  critical 
temperature  gradient  is  exceeded,  the  appropriate  sound 
velocity  will  decrease.  Prom  our  measurement  of  the  deriva¬ 
tives  of  sound  velocity,  we  can  form  the  derivative,  at 
constant  velocity  (P  or  S) 


«n. 


TABLE  III.  Comparison  of  calculated  with  measured  compres- 


I 

I 

I 


sion  for  MgO.11 


Measured 

Compression 

Predicted  from  if 

Murnaghan  Log.  Eq. 

kbars  Measurements 

Polynomial  Eq. 

P 

kbars 

Meas . 

v/v0 

Calc . 

v/v0 

Error,  in 

Percent 

Calc . 

v/vD 

Error,  in 

Percent 

25 

.987 

.986 

-0.10 

.986 

-0.10 

50 

•  97'i 

.972 

-0.21 

.972 

-0.21 

75 

•  963 

.960 

-0.31 

.960 

-0.31 

100 

.951 

.948 

-0.32 

.948 

-0.32 

150 

.930 

.927 

-0.32 

.927 

-0.32 

200 

.910 

.908 

-0.22 

.909 

-0.11 

250 

.893 

.890 

-0.34 

.894 

+0.11 

300 

.877 

.874 

-0.32 

.882 

+0.57 

350 

.862 

.860 

-0.23 

.877 

+1.74 

_ 
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TABLE 

IV.  Comparison  of 

sion  for  Ai^O^ 

calculated  with  measured 

1  2 

t 

compres- 

Measured 

Compression 

PrerH  r't-.Pfi  frnm  il 

Murnaghan  Log.  Eq. 

kbars  Measurements 

Polynomial  Eq. 

P 

Meas. 

Calc . 

Error,  In 

Calc  . 

Error,  in 

kbars 

V/V0 

VA0 

Percent 

v/v0 

Percent 

63 

0.98 

.976 

-0.4 

•  97  6 

-0.4 

128 

0.96 

•  954 

-0,6 

.954 

-0.6 

192 

0.94 

.935 

-0.5 

.935 

-0.5 

256 

0.92 

.918 

-0.2 

.919 

-0.1 

288 

0.91 

.910 

0 

.911 

+0.1 

j^4 

0.905 

.906 

+0.1 

.907 

+0.2 

TABLE  IV.  Comparison  of  calculated  with  measured  compres¬ 
sion  for  * 1 2 


Measured 

Compression 

Predicted  from  4 

Murnaghan  Log.  Eq. 

kbars  Measurements 

Polynomial  Eq. 

P 

kbars 

Meas . 

v/v0 

Calc . 

v/v0 

Error,  in 

Percent 

Calc . 

v/v0 

Error,  in 

Percent 

63 

0.98 

.976 

.=r 

O 

1 

.976 

-0.4 

128 

0.96 

•  954 

-0.6 

.954 

i 

o 

O', 

192 

0.94 

•  935 

-0.5 

.935 

-0.5 

256 

0.92 

.918 

-0.2 

.919 

-0.1 

288 

0.91 

.910 

0 

.911 

+0.1 

304 

0.905 

.906 

+0.1 

.907 

+0.2 

L 
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(20) 


which  yields  the  temperature  depth  relation  for  these 
materials.  For  the  slope  of  velocity  with  depth  to  change 
sign  and,  in  particular  to  become  negative  in  the  upper 
mantle,  the  earth  must  have  properties  such  that 


(d'T/dP ) 


EARTH 


>  (3T/3P) 


V  (MATERIAL) 


(21) 


in  the  upper  mantle,  and  reverses  with  depth. 

From  tables  I  and  II  we  have  the  following  velocity 
derivatives  for  MgO  and  AlpO^ 


MgO 

A12°3 

-78 . 5°C 

34 . 6°C 

-78 . 5°C  25°C 

( 3vp/3P )>p 

7.57* **10“3 

7.711xlO"3 

5.075xl0~3  5.175x10-3 

# 

( 3Vp/3T)p 

-4.3  xlO-4 

-5.0  xio"4 

-3.0  xio"4  -3.7  *10”4 

## 

(3VS/3DT 

4.07*10“3 

5.351*10"3 

2 . 107x10~3  2 . 207x10” 3 

# 

(3Vg/3T)p 

-3.6  xicT4 

-4.8  xi(T4 

-2.2  xio“4  -3.1  xic“4 

*# 

*  km/sec/kbar 

**  km/sec/°C 


and  from  (20)  and  a  value  of  3.3  km/kbar  for  the  pressure- 
depth  relation  in  the  upper  mantle,  we  have  the  results 
(where  Z  is  depth) 
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MgO 

-78 . 5°C 

34.6°C 

A1 

-78.5°C 

2°3 

25°C 

(9T/3Z)Vp 

5.33 

4.67 

5.12 

4.23 

°C/km 

( 3T/3Z ) y 

S 

3.42 

2.75 

2.90 

2.15 

°C/km 

We  estimate  3(3T/3Z)/3T  to  be  -5 .9x10“3/km  for  constant 
v  and  v  for  MgO,  and  -8.6  (constant  v  )  and  -7. 2*10“ 3 

^  S  P 

(constant  vs)for  AljO^.  These  values  probably  are  less  at 
higher  temperatures,  but  do  not  alter  the  argument,  as  at 
higher  temperatures  the  critical  temperature  gradient  would 
be  less  than  at  25°C. 

Values  for  (dT/dZ)EAriT-  which  have  been  deduced  pre¬ 
viously  are  6.6°C/km,i3  r°C/km,14  and  7 . 46°C/km. 1 5  It  is 
apparent  that  these  two  materials,  which  are  constituents 
of  the  earth's  crust,  meet  the  condition  necessary  for  a 
negative  velocity  slope.  Figures  2A  and  2B  illustrate  this 
condition.  In  figures  2A  and  2B  the  temperature-pressure 
profiles  for  two  geothermal  models  are  superposed  on  a  cross¬ 
plot  of  temperature  versus  pressure  at  constant  velocity, 
derived  from  the  data  above.  Figure  3  is  the  result  showing 
how  the  sound  velocities  for  P  and  S  would  vary  with  depth 
for  MgO  and  A^O^  if  the  temperature-pressure  profile  for 
the  earth  is  imposed  on  either  material.  The  relation 
between  the  material  property  and  the  geothermal  gradient  in 


TEMPERATURE  (°c) 


PRESSURE  (KILOBARS) 


’FIG  2A:  T-P  PLOT  AT  CONSTANT  SOUND  VELOCITY  FOR  MgO 
WITH  GEOTHERMIC  GRADIENTS  SUPERPOSED 


TEMPERATURE  (°C) 


PRESSURE  (K!!  OBARS) 


FIG  2B:  T-P  PLOT  AT  CONSTANT  SOUND  VELOCITY  FOR  Al203 
WITH  GEOTHERMIC  GRADIENTS  SUPERPOSED 


NED  FOR  MgO  8 
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producing  a  velocity  minimum  is  made  strikingly  evident.  It 
is  apparent  that  for  a  given  geothermal  gradient,  the  material 
property  of  (3T/3P)y  determines  the  depth  of  the  velocity 
minimum.  Thus  a  known  variation  in  depth  of  the  minimum 
velocity  layer  may  be  attributed  to  variation  in  the  geother¬ 
mal  gradient,  or  material. 

On  the  basis  of  the  above,  it  would  appear  that 

1.  The  observed  velocity  minimum  in  the  mantle 
can  exist  as  a  consequence  of  the  individual  material  pro¬ 
perty,  and  one  need  not  assume,  a  priori,  a  need  for  a  non- 
homogeneous  mantle  to  explain  this. 

2.  Variation  of  the  depth  at  which  the  velocity 
minimum  occurs  may  be  attributed  to  a  variation  in  the  mater¬ 
ial  property. 

3.  Variation  of  velocity  at  constant  depth  could 
be  due  to  composition  difference,  but  as  lateral  variation 
of  the  geothermal  gradient  could  also  produce  this  effect, 
it  may  be  attributable  to  both  factors. 

4.  As  the  critical  velocity  gradients  for  P  are 
greater  than  for  5,  shallow  P  and  pronounced  S  velocity 
minima  result  from  the  behavior  of  the  materials. 

IV.  RECOMMENDATION  FOR  FUTURE  WORK 

The  results  of  the  present  program  have  borne  out  the 
point  of  view,  expressed  in  the  proposal,  that  it  is  higoly 
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profitable  to  perform  experiments  with  precision  enough  to 
define  the  pressure  and  temperature  derivatives  of  the 
velocity  of  mono-mineral  compacts.  Because  of  the  complexity 
of  the  earth,  it  is  important  to  be  able  to  obtain  sufficient 
data,  so  that  the  common  denominators  may  be  abstracted;  and 
from  these  denominators  obtain  a  broader,  and  yet  more 
detailed  understanding  of  the  deep  earth.  We  believe  that 
it  is  Important  to  extend  these  measurements  to  other 
minerals,  either  as  polyeryccalline  compacts  or,  where  these 
are  not  available,  to  single  crystals.  It  is  important, 
for  example,  to  determine  dB/dP  for  other  minerals  and  see 
if  it  is  near  4;  and  if  d2B/dP2  is  zero.  Also  of  importance 
are  the  ways  in  which  the  derivatives  of  velocity  differ  from 
material  to  material,  and  the  effects  of  these  differences 
on  velocities  with  depth.  These  considerations  raise  the 
following  questions: 

1.  Are  the  derivatives  of  such  numerical  values 
that  the  velocities  of  different  minerals  tend  to  converge 
or  diverge  with  increasing  pressure  and  temperature? 

2.  Do  all  oxides  and  silicates  have  critical 
thermal  gradients  similar  to  MgO  and  A^O^? 

These  are  significant  and  important  questions,  which 
bear  directly  on  an  understanding  of  the  mantle.  It  is 
desirable  to  employ  materials  which  are  thought  representa¬ 
tive  of  the  mantle.  However,  this  should  not  be  an  absolute 
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criteria,  f' or  chiocs^Ln^  specimens*  A s  17101*0  information  is 
gained  about  many  materials,  it  will  become  possible  to 
predict  the  behavior  of  any  particular  mineral. 
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APPENDIX  A 
ABSTRACTS  OF  PAPERS 

The  following  abstracts  are  from  the  manuscipts  which 
have  resulted  from  our  effort  in  the  fullfillment  of  our 
contract . 

1.  "Conditions  for  a  density  minimum  in  the  upper  man 
tie",  0.  L.  Anderson,  J .  Geophys .  Res . ,  '£0  [6],  1457-1461, 
1965. 

Abstract  A  number  of  recent  articles  have  indic¬ 
ated  that  there  is  a  density  minimum  in  the  upper 
mantle.  An  exarole  is  the  article  by  .lark  and  Ring- 
wood,  where  the  density  is  based  upon  petrology.  In 
this  paper,  the  conditions  for  a  density  minimum  are 
established  using  the  Mie-Guneisen  equation  of  state. 

It  is  shown  that  for  a  homogeneous  mantle  with  a 
positive  value  of  thermal  expansivity,  the  sign  of 
the  density  gradient  is  the  same  as  the  sign  of  the 
S  velocity  gradient,  that  is,  d(ln  vs)/d(lnp)>0. 

The  conditions  for  a  densi '  . "  minimum  in  a  nonhomo- 
geneous  mantle  are  also  discussed. 

2.  "The  relation  between  refractive  index  and  density 
of  minerals  related  to  the  earth's  mantle",  0.  L.  Anderson 
and  E.  Schreiber,  J .  Geophys .  Res . ,  IQ  [6],  l463-l471>  1965 

Abstract .  It  is  known  that  the  density  of  many  min¬ 
erals  is  related  to  the  (average)  index  of  refraction 


by  a  linear  law  called  the  Gladstone-Dale  law.  It 
is  shown  that  this  law  is  generally  applicable  only 
to  minerals  whose  anion  is  oxygen  and  whose  mean 
molecular  weight  is  close  to  21.  Another  relation¬ 
ship,  called  Crude' s  law,  which  is  deduced  from 
classical  dielectric  theory,  fits  the  data  just  as 
well  as  the  empirical  linear  law.  The  correlation 
between  density  and  index  includes  the  minerals 
arising  from  various  combinations  of  Si02,  MgO, 

AlgOp  Na20,  and  KgO.  An  implication  is  that  den¬ 
sity  controls  certain  physical  properties  of  oxides, 
independently  of  composition  or  crystal  class. 

This  is  analogous  to  the  rule  found  by  Birch  which 
relates  sound  velocity  to  density  and  holds  for  the 
same  groups  of  minerals. 

3.  "Lattice  d  ynamics  in  geophysics",  0.  L.  Anderson,  Tran., . 

N.  Y.  Acad.  Sci.  ,  Ser.  II,  27  [3],  298-308,  1965. 

Abstract .  An  arbitrary  but  useful  way  to  classify 
solid  physical  properties  is  into  the  classes  of 
lattice  properties  or  defect  properties.  In  lattice 
dynamics,  the  word  lattice  indicates  lattice  pro¬ 
perties  and  the  word  dynamics  indicates  the  impor¬ 
tance  of  thermal  vibrations  of  atoms  to  those 
properties.  Examples  of  lattice  dynamics  properties 
are:  infrared  reflection,  specific  heat,  thermal 
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expansivity,  sound  velocity,  melting  temperature, 
and  thermal  diffuse  x-ray  scattering.  These  pro¬ 
perties  are  defined  in  terms  of  the  set  of  normal 
vibrational  modes,  calculated  from  a  model  composed 
of  discrete  masses  connected  by  springs  in  a  large 
array.  This  information  is  contained  in  a  funda¬ 
mental  diagram,  associated  with  the  particular 
model,  called  the  frequency  vs.  wave  number  dia¬ 
gram.  ( v-k) . 

The  importance  of  lattice  dynamics,  as  far  as 
geophysics  is  concerned,  is  that  all  lattice  pro¬ 
perties  are  functions  of  the  atomic  mass,  spacing, 
and  spring  constant  of  the  model.  Therefore,  any 
lattice  property  is  derivable  from  other  lattice 
properties.  In  particular  it  is  possible  to  esti¬ 
mate  several  important  lattice  properties  from  the 
sound  velocity  and  density.  To  geophysicists  the 
field  of  lattice  dynamics  provides  a  set  of  prin¬ 
ciples  by  which  thermal,  optical,  and  thermodynamic 
properties  can  be  estimated  from  seismic  data. 

4.  ’'The  Bulk  Modulus  Volume  Relationship  for  Oxide  Compounds 
and  Related  Geophysical  Problems",  0.  L.  Anderson  and  J.  Nafe, 
J.  Geophys.  Res.,  70  [16],  3951-3963,  1965 
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Abstract .  The  relationship  between  the  sound  velo- 
city  and  density  in  various  oxide  compounds  at- atmo¬ 
spheric  pressure  is  relevant  to  problems  of  «he 
earth's  interior.  Here,  data  on  elastic  constants 
of  various  compounds  are  collected  and  analyzed 
It  is  shown  that  the  bulk  modulus-volume  per  ion 
pair  relationship  for  oxide  compounds  differs  in  a 
remarkable  degree  from  that  found  for  alkali  halides, 
fluorides,  selenides,  sulfides,  and  covalent  com¬ 
pounds  . 

It  Is  shown  that  a  change  of  volume  has  the 
same  effect  on  the  bulk  modulus  of  oxide  compounds, 
whether  the  volume  change  is  produced  by  pressure, 
compositional  variation,  phase  changes,  temperature, 
or  porosity.  It  thus  appears  that  volume  is  the 
primary  variable  affecting  the  elastic  moduli  of 
oxide  compounds,  and  all  other  variables  affect  the 
moduli  only  insofar  as  they  affect  the  volume  itself. 

5.  "An  approximate  method  of  estimating  shear  velocity  from 
specific  heat",  0.  L.  Anderson,  J.  G.  R. ,  70  [18],  ^726-^728, 
1965. 

Abstract .  The  shear  velocity  may  be  estimated  from 
low  temperature  specific  heats  because  at  low  tem¬ 
peratures  the  contribution  to  the  specific  heat  is 
due  to  acoustic  vibrations.  The  low  temperature 
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specific  heat  is  characterized  by  the  "Debye"  tem¬ 
perature  and  the  shear  velocity  may  be  estimated 
from  the  Debye  temperature.  The  source  of  discre¬ 
pancies  and  the  general  implications  are  discussed. 

6.  "The  pressure  derivatives  of  the  sound  velocities  of  poly¬ 
crystalline  magnesia",  0.  L.  Anderson  and  E.  Schreiber,  J.  G.  R. 
70..  [20],  52*11-5248,  1965. 

Abstract .  A  poly  crystalline  sample  of  MgO  was  ob¬ 
tained  of  gem  quality.  The  sound  velocities  and 
the  pressure  derivatives  of  the  sound  velocities 
were  measured  by  the  "phase  comparison"technique . 

The  pressure  derivatives  found  from  experiments  up 
to  4  kilobars  are: 

dvs 

—  -  4.351  x  10  3  km/sec/kbar 

dp 

% 

—  =  7.711  x  10  3  km/sec/kbar 

dp 

Arguments  are  presented  which  indicate  that  these 
pressure  derivatives  hold  up  to  at  least  100  kilo- 
bars  . 

These  values  yield  a  vanishingly  small  value 
of  the  pressure  derivative  of  Poisson’s  ratio. 

The  Gruneisen  constants  of  the  shear  and  long¬ 
itudinal  modes  are  calculated,  from  which  the 


acoustic  Gruneisen  constant  Is  estimated  to  be  1.60. 

This  agrees  very  well  with  the  Gruneisen  constant 
obtained  from  thermal  properties. 

7.  ■  "Two  methods  for  estimating  compression  and  sound  velo¬ 
city  at  very  high  pressures",  0.  L.  Anderson,  Proc .  Nat.  Acad. 
Scl. ,  51  [3],  (367-673,  1965  . 

Abstract .  Precision  sound  velocity  measurements  at 
relatively  low  pressure  (less  than  10  kilobars)  can 
be  used  to  estimate  sound  velocity  and  volume  com¬ 
pression  at  very  high  pressures.  If  the  parameters 
of  a  continuous  analytic  function  are  sufficiently 
accurately  determined  at  low  values  of  the  indepen¬ 
dent  variable,  the  function  may  be  evaluated  at 
high  values  of  the  independent  variable.  Consequently, 
only  materials  undergoing  compression  without  phase 
change  are  treated. 

What  is  new  is  the  presentation  of  two  funct¬ 
ions  suitable  for  representing  volume  dependence  upon 
pressure  (compression  curves)  and  the  evaluation  of 
the  parameters  in  these  functions  from  acoustic  data 
taken  at  low  pressures  but  with  high  precision. 

The  parameters  are  the  bulk  modulus  and  its  higher 
derivatives . 

8.  "The  pressure  derivatives  of  the  sound  velocity  of 
polycrystalline  alumina",  E.  Schreiber  and  0.  L.  Anderson, 

J .  Am.  Ceram.  Soc .  ,  in  press. 
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Abstract .  The  sound  velocities  and  the  pressure 
derivatives  of  the  sound  velocities  were  measured 
on  a  small  sample  of  alumina  (Lucalox)  by  the 
method  of  "pulse  superposition".  The  pressure 
derivatives  found  from  experiments  up  to  b  kilobars 
at  25°C  are: 

dvs 

«  2.207  x  10  3  km/sec/kbar  (Shear  Wave) 
dP 

dvP 

—  =  5*175  x  io  3  km/sec/kbar  (Long.  Wave) 
dP 

These  values  lead  to  the  variation  of  Poisson’s 
ratio  with  pressure  of  da/dP  =  1.02  *  10“4/kilo- 
bar;  and  to  the  variation  of  (isothermal)  bulk 
modulus  with  pressure  of  dB^/dP  *  3*99* 

Using  the  above  values  the  computed  compress¬ 
ion  is  determined  up  to  *100  kilobars,  and  compares 
well  with  shock  wave  measurements.  Consequently, 
the  measured  pressure  derivatives  may  hold  up  to 
several  hundred  kilobars. 

9.  "Seismic  parameter  Computation  at  very  hig  .  pressure 
from  laboratory  data:,  0.  L.  Anderson,  Bull.  Seism.  Soc.  Am., 
in  press. 

Abstract .  By  using  the  accuracy  inherent  in  ultra¬ 


sonic  velocity  measurements  taken  at  pressures  less 
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than  10  kb,  the  seismic  parameter  $  *  v*  -  (V3)  v§ 
can  oe  computed  at  very  high  pressures.  The  equa¬ 
tion  used  requires  the  assumption  that  the  second 
derivative  with  respect  to  pressure  of  the  bulk 
modulus  be  negligible  at  all  pressures  considered. 

This  assumption  is  checked  by  computing  the  compres¬ 
sion  (V/V  )  in  the  pressure  range  by  equations  of 
state  using  the  assumption,  and  comparing  the  res- 
sulting  values  with  measured  compression.  Illust¬ 
rations  are  given  for  MgO  and  A^O^. 

10.  "Temperature  dependence  of  the  velocity  derivatives  of 
periclase",  E.  Schreiber  and  0.  L.  Anderson,  submitted  for 
publication. 

Abstract .  The  temperature  dependence  of  the  sound 
velocity  in  polycrystalline  MgO  has  been  determined 
from  +80  to  -80°C  and  to  2  kbars  at  -78.5°C. 

These  results  are  compared  to  previous  measurements 
to  4  kbars  at  3^.6°C.  From  these  measurements  the 
critical  temperature  gradient  (3T/sP)v  for  MgO  was 
determined,  and  these  values  applied  to  the  earth. 
Because  (3T/3P)v  was  found  to  be  greater  for  P 
waves  than  for  S  waves  in  MgO,  and  both  are  small 
compared  to  (3T/9P)V  for  the  earth,  it  is  possible 
to  explain  the  existence  of  a  low-velocity  layer 
in  a  homogeneous  earth. 
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APPENDIX  B 

MEASUREMENT  OF  SOUND  VELOCITY  BY  MEANS  OF 
ULTRASONIC  INTERFEROMETRY 

A .  Phase  Comparison 

The  technique  involves  many  of  the  optical  principles 
of  multiple-beam  interferometry.  The  ultrasonic  technique, 
called  "phase  comparison"  method,  was  originated  by  McSkimin1 
and  is  described  below. 

Figure  1A  illustrates  the  electronic  arrangement.  A 
stable,  variable  frequency  oscillator  is  employed  as  the 
RF  carrier  generator.  This  frequency  is  accurately  measured 
by  the  frequency  counter.  The  output  of  this  oscillator  is 
fed  to  the  gated  harmonic  generator.  The  purpose  of  this 
component  is  to  produce  series  of  gated  RF  pulses  at  the 
desired  frequency  (in  this  case,  60  mc/sec).  The  frequency 
and  width  of  these  gated  pulses  are  controlled  by  the  General 
Radio  1217-B  unit  pulse  generator.  The  output  of  the  harmonic 
generator  is  applied  directly  to  a  quartz  transducer  (x  or  y 
cut  for  longitudinal  or  shear  modes)  which  serves  as  a  trans¬ 
mitter  of  the  generated  pulse  and  receiver  of  the  echoes  in 
the  time  slot  between  the  generated  pul  s.  The  beat  fre¬ 
quency  oscillator  is  used  for  heterodyning  to  produce  a  20  me 
I.F.  thus  permitting  the  use  of  a  single  high  gain  I.F. 
amplifier.  The  output  is  also  tied  to  an  attenuator  and  then 
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to  the  high  gain  I.P.  amplifier.  The  output,  of  this  ampli¬ 
fier  is  displayed  on  a  Tektronix  535A  oscilloscope. 

To  measure  the  absolute  velocity  at  ambient  conditions 
a  buffer  rod  of  either  quartz  or  vitreous  silica  is  employed. 
A  quartz  transducer  is  bonded  to  one  end  of  the  buffer  rod 
and  performs  as  the  driver.  The  sample  is  bonded  to  the 
other  end  of  the  rod  with  a  suitable  material  (Dow-Corning 
resin  blend  V-9  for  example).  Part  of  the  energy  trans¬ 
mitted  by  the  transducer  is  reflected  at  the  opposite  end 
of  the  buffer  rod,  and  part  is  coupled  into  the  specimen, 
via  the  bond,  to  be  reflected  back  from  the  free  end  of  the 
specimen.  The  output  seen  on  the  oscilloscope  consists  of 
the  pulse  applied  to  the  transducer  from  the  harmonic  gener¬ 
ator,  the  return  echoes  from  the  end  of  the  buffer  rod;  and 
between  the  buffer  rod  echoes,  the  specimen  echoes.  The 
energy  coupled  from  the  buffer  rod  is  now  the  applied  pulse 
to  the  specimen.  By  broadening  the  pulse  width  so  that  the 
energy  returning  from  the  specimen  echoes  is  overlapped  by 
the  applied  pulse  from  the  buffer  rod,  the  60  me  carrier 
is  made  to  phase  interfere.  It  is  in  this  manner  that  the 
accuracy  is  obtained.  A  minima  will  be  obtained  when  the  RP 
frequency  is  such  that  specimen  length  is  an  integral  number 
of  wave  lengths.  Prom  the  condition  for  interference,  the 
velocity  is  given  by1 
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v 


2  i  f 
n  +  9/360 


(1) 


where  v  is  the  velocity  in  the  specimen,  i  is  the  specimen 
length,  f  th j  f quency ,  n  the  number  of  wave  lengths  in 
the  specimen,  and  e  a  small  correction  for  phase  shift  intro¬ 
duced  in  the  seal.  The  latter  may  be  evaluated  from  the 
known  acoustic  impedances  of  the  buffer  rod,  seal,  and  spe¬ 
cimen  as  described  by  McSkimin1.  When  this  correction  is 
applied,  velocity  can  be  determined  to  an  accuracy  of  at 
least  one  part  in  10,000.  The  order  n  is  determined  by 
measuring  a  sequence  of  frequencies  at  which  'resonance* 
occurs  and  by  using  the  following 

f 

—  -  ft  /  ^  \ 


where  fn  is  the  resonant  frequency  for  a  particular  n,  and 
Af  the  difference  in  ffi  and  f^n  + 

To  measure  the  velocity  with  pressure  or  temperature, 
the  transducer  is  bonded  directly  to  the  specimen  and  the 
specimen  echoes  are  overlapped.  The  change  in  frequency  at 
constant  n  is  determined.  The  velocity  at  any  pressure  or 
temperature  may  be  calculated  from 


v  ■  vo<fr/fn0>(l/to> 


(3) 
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where  the  o  subscripts  refer  to  initial  conditions.  The 
ratio  may  be  obtained  from  expansivity  data  for  temper¬ 

ature  runs.  The  situation  is  slightly  more  involved  for 
pressure  runs.  Cook2  has  shown  that  the  value  of  the  ratio 
£/£q  may  be  determined  from  the  relation 


o  _  ,  .  1  +  A  !  d  P 

i  '  1  p  4  a  -  3B 

o  dQ 


(4) 


where  A  is  given  by  (vQfn/fn  )2  for  longitudinal  waves  and 

B  by  the  value  (vQf  /f  )2  for  shear  waves.  The  initial 

n  n0 

density  is  p0  and  A  is  given  by 


where  Cp  is  the  heat  capacity,  0  is  the  coefficient  of 
linear  expansion,  Ba  is  the  adiabatic  bulk  modulus,  and  T 
is  the  absolute  temperature.  The  atmospheric  pressure  values 
of  the  parameters  are  used  in  determining  A,  which  leads  to 
errors  no  more  than  about  0.01$,  and  the  value  of  lQ/i  is 
found  as  a  function  of  pressure.  In  this  manner  the  change 
of  path  length  with  pressure  is  accurately  accounted  for. 


B.  Pulse  Superposition 


Pulse  superposition  differs  from  phase  comparison  in 
that  the  frequency  of  the  applied  pulses  is  crucial  to  the 
measurement  rather  than  the  frequency  of  the  carrier.  This 
technique  is  also  due  to  McSkimin.3  It  has  the  operating 
advantage  of  greater  energy  in  the  return  echoes.  Figure  IB 
is  the  block  diagram  of  the  electronic  arrangement  for 
pulse  superposition  (P.S.P.).  The  critical  component  is  the 
pulsed  oscillator .#  The  C.W.  oscillator  is  of  variable 
frequency  design,  with  a  range  of  8  to  130  mc/sec  permitting 
the  use  of  a  wide  range  of  transducers.  The  C.W.  oscillator 
is  pulsed  by  the  pulse  repetition  oscillator  at  a  pulse 
repetition  frequency  (P.R.F.)  of  1  me  or  less.  The  output 
of  the  Uni-Pulse  consists  of  a  sequence  of  R.F.  pulses  which 
drive  the  X  or  Y  cut  transducers.  The  returning  specimen 
echoes  are  received  by  the  same  transducer.  The  signal  is 
fed  to  a  mixer  stage,  a  calibrated  attenuator,  and  a  high 
gain  I.F.  amplifier-detector.  The  detected  signal  Is 
applied  to  a  Tektronix  535A  oscilloscope,  where  the  envelopes 
of  the  pulses  are  displayed.  The  P.R.F.  is  accurately  deter¬ 
mined  to  six  significant  figures  with  the  frequency  counter. 
This  method  requires  a  very  stable  P.R.F.,  at  least  1  part 
in  107. 

#Uni-Pulse,  mftrd  by  the  7.M.A.  Anderson  Co.,  Pleasant 
Valley,  Conn. 
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In  principle,  the  P.R.F.  is  adjusted  so  that  Its  period 
is  equal  to  some  integral  sub-multiple  of  the  delay  time  in 
the  specimen.  That  is,  the  time  delay  between  applied  pulses 
(period  of  P.R.F.)  is  exactly  equal  to  an  integral  number  of 
round  trips  in  the  specimen.  When  this  condition  is  achieved, 
the  applied  pulses  are  superposed  upon  the  specimen  echoes. 

If  the  Integer  is  one,  every  specimen  echo  will  have  an 
applied  pulse  superposed  upon  it.  If  the  integer  is  two,  the 
applied  pulses  will  be  superposed  upon  every  other  echo,  and 
so  on.  It  is  preferable  to  operate  with  this  integer  equal 
to  unity,  for  then  the  greatest  amount  of  energy  is  being 
impressed  upon  the  specimen.  At  this  condition,  only  the 
applied  pulses  are  visible  on  the  oscilloscope  c' splay.  To 
observe  the  specimen  echoes,  a  "window"  is  produced  in  the 
sequence  of  applied  pulses  by  gating  the  pulsed  oscillator. 
This  is  accomplished  by  applying  a  gating  voltage  from  the 
oscilloscope  to  the  Uni-Pulse.  One  now  views  the  echoes  and 
critically  adjusts  the  P.R.F.  to  maximize  their  amplitudes. 
When  critically  adjusted,  the  applied  pulse  is  superposed 
upon  the  echo,  and  the  C.W.  in  the  pulse  is  phase  adding 
with  the  C.W.  in  the  echo.  The  relation  between  the  measured 
time  delay  (reciprocal  of  the  P.R.F.)  and  the  actual  time 
delay  in  the  specimen  is  given  by3 


(5) 
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In  (1)  t  Is  the  measured  period  of  the  P.R.F.,  at  the 
interference  condition,  6  is  the  true  time  de)ay  in  the 
specimen,  p  is  the  integer  discussed  above,  ?  is  the  phase 
shift  introduced  by  the  seal  (transducer  to  specimen  bond), 
f  is  the  frequency  of  the  C.W.  oscillator,  and  n  is  an  integer 
associated  with  the  phasing  between  the  C.W.  within  the 
applied  pulse  and  within  the  return  echo.  To  understand 
the  meaning  uf  n,  consider  Figure  2. 

A  constructive  interference  will  occur  every  time  a 
C.W.  cycle  in  the  pulse  is  exactly  in  phase  with  a  C.W. 
cycle  in  the  echo.  A  series  of  maxima  may  therefore  be 
observed  as  the  P.R.F.  is  varied.  These  maxima  will  be 
separated  by  differences  in  the  P.R.F.  corresponding  to  the 
period  of  the  C.W.  frequency.  The  maxima  will  be  observed 
for  each  Integral  n  for  n  £  0  £  n.  It  is  possible  to  deter¬ 
mine  the  P.R.F.  for  which  n  =  0,  as  described  by  McSkimin.3 

Choosing  the  conditions  p  =  1  and  n  =  0,  the  true  time 

delay  in  the  specimen  (time  per  round  trip)  is  given  by 

6  =  t  +  s/360f  (6) 

Generally,  r,  is  less  than  1°  for  a  properly  prepared 
seal,  and  the  C.W.  carrier  is  of  the  order  of  107 *  cycles/sec . 
The  correction  is  of  the  order  of  10"10  or  less,  so  that 
even  for  precise  measurements  this  correction  may  frequently 
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FIG  2 :  ILLUSTRATING  THE  MEANING 
OF  THE  INTEGER  N 
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be  ignored.  The  simple  manner  of  dealing  with  these  sources 
of  systematic  errors  is  a  major  advantage  of  pulse  super¬ 
position. 

A  further  advantage  arises  from  the  fact,  that  for  C.W. 
frequencies  near  the  transducer  resonance  frequency,  the 
value  of  c  is  sensitive  to  changes  of  temperature  and 
pressure . 

The  velocity  of  sound  in  the  specimen  is  readily  obtained 
once  the  specimen  length  z  is  known.  The  velocity  is  given 
jy 


As  6  is  generally  known  to  5-6  significant  figures, 
the  greatest  limitation  on  the  accuracy  of  v  is  the  measure¬ 
ment  of  the  specimen  length  in  the  direction  of  sound 
propagation. 

The  measurement  of  the  velocity  with  pressure  or  tem¬ 
perature  is  the  same  as  for  phase  comparison.  The  P.R.P. 

(at  constant  n)  is  measured  and  the  ratio  of  f/f0  (<50/6) 

used  in  (3)  in  place  of  fn/f  .  The  length  ratio  is  also 

o 

determined  in  the  same  manner  as  for  phase  comparison. 
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APPENDIX  C 

ADIABATIC  TO  ISOTHERMAL  CORRECTIONS  TO  THE  SECOND  ORDER 

The  results  of  acoustic  experiments  yields  the  adia¬ 
batic  elastic  constants.  It  Is  desirable  for  numerous 
reasons  to  convert  from  adiabatic  to  Isothermal  conditions. 
In  order  to  preserve  the  inherent  accuracy  of  the  acoutis 
measurements,  an  exact  treatment  Is  desirable.  Overton1 
developed  equations  for  cubic  crystals.  Be3ow  these  are 
presented,  modified,  for  the  isotropic  case. 

The  adiabatic-isothermal  corrections  are 


Cp  -  Cv  =  TVa2/ XrT, 


(1) 


XT  -  xs  =  TVa2/Cp 


(2) 


G 


s 


(3) 


where  a  is  the  volume  coefficient  of  thermal  expansion,  Cv 
and  C  represent  specific  heats,  V  represents  specific 

r 

volume,  x  represents  compressibility,  G  represents  shear 
modulus,  and  the  subscripts  s  and  T  indicate  adiabatic  and 
isothermal  conditions. 


We  have  further  that 


XT 

bt  =  Cv  =  xs 
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where  B  represents  bulk  modulus . 

The  Gruneisen  ratio  is  given  by 


T  csV 

Y  *  Cv  =  Cv*T 


so  that  (1)  and  (2)  also  have  the  forms 


Cv(l  +  Tcy) 


Xs(l  +  Toy) 


Differentiating  (2)  we  have. 


l)  fdxs) 

J  IJ  J  |p 


TVa  2 
C 

P 


.2/12. 

XT  +  a  (  3P 


_i_  (\ ®£pA " 
Cp  l  3  P  J  ip 


We  now  change  the  terms  on  the  right  side  of  (7)  to  measur¬ 
able  quantities.  Using  the  identity. 
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we  have 

j 


Similarly 


Equation  (7)  becomes 


The  equivalent  of  (10)  in  terms  of  bulk  modulus  is,  after 

some  reduction. 

* 
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The  second  term  of  (11)  is  small  and  positive  and  the 
third  term  is  even  smaller  and  ordinarily  negative.  The 
product  Tay  is  usually  small,  of  the  order  of  1  or  2%  at 
ordinary  temperatures. 

The  terms  used  in  obtaining  the  conversion  for  the 
adiabatic  to  isothermal  conditions  for  the  pressure  varia¬ 
tion  of  the  bulk  modulus  are  given  in  Tables  I  through  IV, 
for  both  MgO  and  AljO^  at  the  two  temperature  at  which  the 
isobaric  experiments  were  performed. 
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TABLE  I.  Correction  terms  in  the  calculation  of  (3BT/3P)  . 

_  _  ^ 


i 

i 

PARAMETER 

SYMBOL 

VA'.UE 

UNITS 

i. 

Adiabatic  P. 

5  * 

deriv.  of  B 

OBs/3P)T 

3.916 

% 

Volume 

i 

expansivity 

a 

31.4xl0“6 

i 

Specific  heat 

C 

P 

9.4l*106 

erg/gm/°C 

Adiabatic 

*  - 

compressibility 

xs 

5.82i»xio“13 

cm2/dyne 

I 

Thermal  Gruneisen 

* 

constant 

Y 

1.585 

r 

Isothermal 

I 

compressibility 

5 .913x10“l 3 

cm2/dyne 

Temp,  deriv. 

bulk  modulus 

(3BT/3T)p 

0.291 

kbar/°K 

Temp,  deriv. 

r 

expansivity 

( 3a/3T)p 

IxlO"7 

i- 

Adiabatic- 

| 

isothermal corr , 

1+ayT 

1.0153 

£ 

Coefficient  of 

r 

[A-B],  Eq.  (11) 

ayT(Cv/Cp) 

0.0151 

L 

Value  of  A, 

«• 

Eq.  (11) 

2(aBT)-l(3BT/3T)p 

-10-96 

L 

Value  of  D, 

Eq.  (11) 

a-2(3  x/3T)p 

103 

r 

L 

2nd  term  of 

Eq.  (11) 

0.0M 

” 

1 

3rd  term  of 

«*** 

Eq.  (11) 

-0.022 

Isothermal  P. 
deriv.  of  B 


(3BT/3P)T 


3.94 


TABLE  II.  Corrections  for  computing  (3BT/3P)  of  MgO  at 
-78. 5°C 


PARAMETER 

SYMBOL 

VALUE  UNITS 

*P.  deriv. 
adiabatic,  Bg 

(  9  Bg/  3  P  )ip 

4.00 

*Vol.  expansivity 

av 

22 .08*10“6  °C~1 

:5Specific  heat 

C 

p 

6.37x10s  erg/gm/°C 

#adiabatic 

compressibility 

*S 

5,7823*10"13  cm2/gm 

+Gruneisen 

constant 

Y 

1.67 

■(■Isothermal 

compressibility 

Xiji 

5 • 8238X10"1 3  cmz/gm 

“T.  deriv., 
bulk  modulus 

(3BT/3T)p 

-0.25  kbar/°K 

*T.  deriv. , 
expansivity 

(3av/3T)p 

1 . 3X10*"7 

tAdiabatic- 

isothermal  corr. 

1+ovyT 

1.00717 

vCoeff . :  2nd  term, 
Eq.  (11) 

avyT(Cv/Cp) 

7.11xl0"3 

+2nd  term  of  2nd 
term,  Eq.  (11) 

2(aBT)-i(3ET/3T)p 

-12.9 

r3rd  term  of  3rd 
term,  Eq.  (11) 

av-2(3av/3Dp 

266 

tTotal  2nd  term, 

Eq.  (11) 

o  .034 

tTotal  3rd  term, 

Eq.  (11) 

-0.013 

tlsothermal  P 
deriv.,  bulk  mod. 

(  3B^i/3P  )rp 

4.02 

^Measured 

tCalculated 


TABLE  III.  Corrections  for  computing  (3Brp/3P)  of  Al„0o  at 
25°C.  1  5 


PARAMETER 

SYMBOL 

VALUE 

UNITS 

^Adiabatic  P. 
deriv.,  B_ 

( 3  Eg/  3  P )  ip 

3-98 

*V'ol.  expansivity 

a 

16 . 32*10“6 

°c-1 

#Specific  heat 

°p 

7.831xl06 

erg/gm/°C 

#Adiabatic 

compressibility 

3.967*10~13 

cm2/gm 

tGruneisen  const. 

Y 

1.323 

tlsothermal 

compressibility 

Xip 

3-993x10-13 

cm2/gm 

WT.  deriv. , 
bulk  modulus 

(3BT/3T)p 

0.207 

kbar/°K 

*T.  deriv., 
expansivity 

(3a/3T)p 

1.6l*10“8 

tAdiabatic- 

isothermal  corr. 

1+ayT 

1.006437 

tCoeff.j  2nd  torm 

ayT(C V/Cp) 

6.395x10"3 

t2nd  term  of  2nd 
term,  Eq.  (11) 

2(aBT)-1(3ET/3T)p 

-10.129 

+3rd  term  of  3rd 
term,  Eq.  (11) 

a”2 (3a/3T)p 

60.44 

tTotal  2nd  term, 

Eq.  (11) 

0.012 

tTotal  3rd  term, 

Eq.  (11) 

-0.002 

tlsothermal  P. 
deriv.,  bulk  mod. 

( 3  Bip/  3  P  )  ip 

3.99 

^Measured 

tCalculated 


TABLE  IV.  Corrections  for  computing  (3B-/3P)  of  A1  0  at 
-78.5°C.  t  2  3 

PARAMETER  SYMBOL  VALUE  UNITS 


^Volume 

expansivity 

ct 

9.21xl0“6 

#P.  deriv. , 

adiabatic,  B 
’  s 

(9BS/3P)T 

4.00 

^Specific  heat 

C 

P 

4.8l6xl06  erg/gm/°C 

^Compressibility , 
adiabatic 

xs 

3. 9460x10“ 13  em?-/gm 

tGruneisen 

constant 

Y 

1.226 

tlsothermal 

compressibility 

XT 

3*9547x101 3  cm2/gm 

*T.  deriv., 
bulk  mod. 

OBT/aT)p 

0.1284  kbar/°K 

*T .  deriv., 
expansivity 

(9a/3T)p 

2 .73x10”8 

tAdiabatic- 

isothermal  corr. 

1+ayT 

1.002213 

tCoeff.;  2nd  term 
Eq.  (11) 

ayT(Cv/Cp) 

2.208xl0”3 

+2nd  term  of  2nd 
term,  Eq.  (11) 

2(aBT)_1(3BT/3T)p 

-11.02 

+3rd  term  of  3rd 
term,  Eq.  (11) 

a“2 ( 3a/3T)p 

321.8 

tTotal  2nd  term, 
Eq.  (11) 

0.007 

tTotal  3rd  term, 
Eq.  (11) 

0.006 

tlsothermal  P. 
deriv.  bulk  mod. 

( 3B(j/3P  )fy 

4.01 

-  4H 
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APPENDIX  D 

EXPERIMENTAL  DATA 

The  measurements  involve  the  determination  of  5  time-1), 
a  pseudo-resonant  frequency,  and  the  frequency  ratio  as 
outlined  in  Appendix  B.  Here,  the  data  obtained  in  the 
course  of  the  experiments  and  the  data  employed  in  subse¬ 
quent  computations  are  presented. 

A.  MgO 

The  length  of  the  MgO  specimen  used  was  0. 6692*4  cm  and 
the  bulk  density  3*5803  at  25°C.  The  initial  velocities 
were  vp  *  9*7711  and  vs  =  5*9682  as  determined  at  25°C 
and  1  atm. 

Runs  at  1  atm  and  temperatures  from  80°C  to  -78.5°C 
(the  sublimation  temperature  of  dry  ice)  were  performed  and 
the  velocities  determined  as  a  function  of  temperature.  The 
data  for  these  runs  are  given  in  Table  I.  The  value  of  f  0 / £ 
was  obtained  from  the  thermal  expansion  data  of  MgO.  This 
has  been  evaluated  by  White  and  Anderson,1  based  on  the 
work  of  Ganesan,2  Sharma,3  and  Skinner.4 

Variation  of  the  velocities  with  pressure  at  constant 
temperature  were  performed  at  +3*4. 6°C  and  -78.5°C.  The 
vessel  was  thermostated  with  a  heater  at  +3*4. 6°C  and  with 
dry  ice  at  -78.5°C.  The  data  and  velocities  obtained  are 


TABLE  I, 

Longitudinal  and  Shear  Velocities 
perature  at  1  atm  for  MgO. 

with  Tem- 

^  Temp. 

I  °C 

z0/i 

Frequency  Ratio 

Long.  Shear 

Velocity 

Long. 

(km/sec) 

Shear 

I  -8t> 

0.999046 

1.005744 

1.008665 

9.8186 

6.0145 

-78.5 

0.999057 

1.005669 

1.008591 

9. 3179 

6.0139 

1  -70 

0.999120 

1.005232 

1.008083 

9.8143 

6.0109 

-50 

0.999200 

1.004?15 

1.007340 

9.8100 

6.0069 

4  -30 

0.995282 

1.004193 

1.006520 

9.8057 

6.0027 

-40 

0.999368 

I.OC3663 

1.005666 

9.8014 

5.9982 

-30 

0.999457 

1.003115 

1.004835 

9.7969 

5.9938 

1  -20 

0.999549 

1.002555 

1.003987 

9.7924 

5.9893 

-10 

0.999644 

1.001981 

1.003130 

9.7877 

5.9847 

0 

0.999742 

1.001401 

1.002251 

9.7830 

5.9801 

a** 

10 

0.999843 

1.000819 

1.001365 

9.7783 

5.9754 

20 

0.999947 

1.000234 

1.000460 

9.7736 

5.9706 

1  25 

i.000000 

1.000000 

1.000000 

9.7711 

5.9682 

I  30 

1.000053 

0,999691 

0.999546 

9.7686 

5. 965? 

40 

1.000159 

0.999072 

0.998638 

9.7636 

5.9610 

50 

1.000268 

0.998453 

0.997730 

9.7586 

5.9562 

60 

1.000378 

0.997834 

0.996822 

9.7536 

5.9515 

1  70 

1.000491 

0.997215 

0.995914 

9.7487 

5.9467 

'f  80 

1.000604 

C. 996597 

0. 995006 

9 . 7'*37 

5.9420 
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listed  in  Taole  II.  Table  III  lists  the  variation  of  bulk 
modulus  with  pressure  computed  from  these  data  at  the  two 
temperatures . 

s.  ai2o3 

The  experiments  on  parallel  those  for  MgO.  The 

length  of  the  A1203  specimen  used  was  1.2340  cm  and  the  bulk 
density  3-972  g/cm2  at  25°C.  The  initial  velocities  were 
determined  to  be  10.845  km/sec  for  the  longitudinal  velocity 
and  6.3730  km/sec  for  the  shear  velocity. 

The  variation  of  velocity  with  temperature  as  determined 
from  25°C  to  -78.50C  is  given  in  Table  IV.  The  thermal 
expansion  data  of  Wachtman5  was  employed. 

Table  V  lists  the  velocity  data  obtained  for  A120^  as  a 
function  of  pressure  at  constant  temperature.  Table  VI 
lists  the  bulk  modulus  as  a  function  of  pressure  at  the  two 
temperatures . 


TABLE  II.  Longitudinal  and  shear  velocities  with  pressure 
at  -78.5°C  and  34.6°C. 


Pressure 

baj 

V* 

Frequency  Ratio 
Long.  Shear 

Velocity 

Long. 

(km/sec) 

Shear 

34.6°C 

1 

1.0000000 

1.000000 

1.000000 

9.7662 

5.5635 

500 

1.0000973 

1.000492 

1.000463 

9.7700 

5.9657 

1000 

1.0001946 

1.000985 

1.000927 

9.7739 

5.9679 

1500 

1.0002919 

1.001478 

1.001390 

9.7778 

5.9700 

2000 

1.0003892 

1.001970 

1.001853 

9.7816 

5.9722 

2500 

1.0004865 

1.002463 

1.002317 

9-7855 

5.9744 

3000 

1.0005838 

1.002955 

1.002780 

9.7893 

5.9766 

3500 

1.0006811 

1.003448 

1.003243 

9-7932 

5.9787 

4000 

1.0007784 

1.003940 

1.003706 

9.7970 

5.9809 

-78.5°C 

1 

t 

1 

1.0000000 

1.000000 

1.000000 

9.8179 

6.0139 

500 

1.0001096 

1.000496 

1.000451 

9.8216 

6.0158 

1000 

1.0002192 

1.000992 

1.000903 

9.8254 

6.0179 

1500 

1.0003288 

1.001488 

1.001355 

9.8292 

6.0199 

20u^ 

1.0004383 

1.001984 

1.003806 

9.8330 

6.0220 

TABLE  III.  Variation  of  adiabatic  bulk  modulus  with  pressure 
at  34.6°C  and  -78.5°C  for  MgO. 

Pressure  Density  Bulk  Modulus  (kbars) 


(bars ) 

34.6°C 

-78 . 5°C 

34.6°C 

-78 , 5°C 

1 

3.5792 

3-5903 

1716.7 

1729.4 

500 

3-5802 

3.5914 

1719.7 

1731.4 

1000 

3.5813 

3.5926 

1720.5 

1733-5 

1500 

3-5823 

3.5938 

1722.6 

1735.6 

2000 

3.5834 

3.5939 

1724.5 

1737.6 

2500 

3.5844 

1726.4 

3000 

3.5855 

1728.4 

3500 

3.5865 

1730.4 

4000 

3.5875 

1732.3 

1 

m 


TABLE  IV.  Longitudinal  and  shear  velocity  as  function  of 
temperature  at  1  atm. 


Temp. 

l/% 

Frequency  Ratio 

Velocity 

(km/sec) 

°C 

0 

Long. 

Shear 

Long. 

Shear 

I 

1 

1- 

25 

1.000000 

1.000000 

1.000000 

10.845 

6.3730 

20 

0.999985 

1.000186 

1.000272 

10.847 

6.3746 

Sr 

10 

0.999929 

1.000567 

1.000807 

10.850 

6.3777 

a-3' 

| 

3 

* 

0 

0.999879 

1.000946 

1.001332 

10.854 

6.3807 

j-j. 

-10 

0.999830 

1.001320 

1.001847 

10.857 

6.3837 

f 

1- 

-20 

0.999782 

1.0016'  ; 

1.002362 

10.861 

6.3866 

r 

-30 

0.999736 

1.0021  ,0 

1.002872 

10.864 

6.3896 

-HO 

0.999697 

1.002423 

1.003367 

10.868 

6.3925 

r 

5 

m--- 

-50 

0.999657 

1.002771 

1.003857 

10.871 

6 . 395*1 

-60 

0.999619 

1.003097 

1.004327 

10.874 

6.3981 

I 

-70 

0.999583 

1.003418 

1.004749 

10.877 

6.4006 

j 

i 

-78.5 

0.999554 

1.003680 

1.005080 

10.8799 

6.40251 

I. 

r 

-80 

0.999549 

1.003722 

1.005138 

10.880 

6.4028 

i 


r 

c 

p 


1 

i= 

TABLE  V. 

Longitudinal 

and  shear 

velocity  as 

function 

of 

5 

i 

i 

pressure  at 

25°C  and  - 

78 . 5°C . 

Pressure 

V* 

Frequency  Ratio 

Velocity 

(km/sec ) 

ii 

(bars ) 

Long. 

Shear 

Long. 

Shear 

i 

i 

i. 

2fj°C 

1 

1.0000000 

1.000000 

1.000000 

10.845 

6.3730 

1 

* 

500 

1.0000658 

1.000305 

1.000239 

10.848 

6.3741 

r 

1000 

1.00013?5 

1.000611 

1.000479 

10.850 

6.3752 

i 

1500 

1.0001982 

1.000917 

1.000719 

10.853 

6.376? 

r 

i 

& 

2000 

1.0002648 

1.001223 

1.000958 

10.855 

6  ?774 

r 

2500 

1.0003304 

1.001528 

1.001198 

10.858 

6.3785 

l 

3000 

1.0003969 

1.001834 

1.001438 

10.861 

6.3796 

| 

3500 

1.0004624 

1.002140 

1.001677 

10.863 

6.3807 

1 

4000 

1.0005300 

1.002446 

1.001917 

10.866 

6.3818 

r 

1 

-78 . 5°C 

r 

i 

1 

1.0000000 

1.000000 

1.000000 

10.8799 

6.40251 

500 

i. 0000656 

1.000300 

1.000230 

10.883 

6.4035 

r 

1000 

1,0001317 

1.000600 

1.000461 

IO.885 

6.4046 

1500 

1.0001973 

1.000900 

1.000691 

10.688 

6.4057 

r 

2000 

1.0002632 

1.001200 

1.000922 

10.890 

6.4067 

r 

2500 

1.0003287 

1.001499 

1.001152 

10.893 

6.4078 

l 

3000 

1.0003947 

1.001799 

1.001383 

10.895 

6.4088 

r 

3500 

1.0004601 

1.002099 

1.001613 

10.898 

6.4099 

L, 

5 

4000 

1.0005254 

1.002399 

1.001844 

10.900 

6.4109 

r 

L. 


TABLE  VI.  Variation  of  Adiabatic  Bulk  Modulus  with  pressure 
at  25°C  and  -78.5°C  for  AlgOg. 


Pressure 
(bars ) 

Density 

25°C 

(gms/cm3) 

-78 . 5°C 

Bulk  Modulus 
25°C 

(kbars ) 
-78.5°C 

1 

3.9720* 

3-9773 

2520.6 

2534.2 

500 

3-9728 

3.9781 

2522.5 

2535.8 

1000 

3.9736 

3.9789 

2524.5 

2538.1 

1500 

3.9744 

3-9797 

2526.5 

2540.1 

2000 

3.9752 

3.9804 

2528.5 

2542.1 

2500 

3-9759 

3.9812 

2530.5 

2544.1 

3000 

3.9767 

3.9820 

2532.5 

2546.1 

3500 

3.9775 

3.9828 

2534.5 

2548.1 

4000 

3.9783 

3.9836 

2536.5 

2550.1 

*Value  of  p 

=  3.972  at 

25°C  and  1 

atm  is  assumed  to 

be 

exact  for  the  purpose  of  this  calculation. 


W  O 
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pressure  derivative  B0',  which  predicts  compression  at  extremely  hig  i 
pressures.  Using  the  derived  expression  Jtn(VD/V)  =  [  1/B0 '  £n(P/B0)+l] 
and  the  values  of  B0  and  B0’  for  MgO  and  AI2O3  reported  here,  the 
agreement  between  the  calculated  compression  and  measured  compression 
reported  in  the  literature  was  within  0.5#  to  350  kilobars.  ; 
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